Generating Function in Quantum Mechanics: An Application to Counting Problems 
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In this paper we present a generating function approach to two counting problems in elementary 
quantum mechanics. The first is to find the total ways of distributing N identical particles among 
m different states. The second is to find the degeneracies of energy levels in a quantum system with 
multiple degrees of freedom. Our approach provides an alternative to the methods in textbooks. 
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I. INTRODUCTION 

Generating function (GF) is a terminology frequently 
used in combinatorial mathematics as a bridge between 
discrete mathematics and continuous analysis. It func- 
tions like a clothesline on which we hang up a sequence 
of numbers for display 1 . A GF may have different aliases 
in different contexts, e.g. partition function in statistical 
physics^ and Z-transform in signal analysis^. It can be 
used to find certain characteristics of a sequence like ex- 
pectation value or standard deviation^. Sometimes that 
will be much more efficient than a direct evaluation. 

There are many ways to construct a GF for a sequence, 
the simplest of which is using power series. Given {a n } 
(n > 0), we may construct its GF as 



(1) 



n>0 



Note that t can take a continous value and F(t) is ana- 
lytic with respect to t, which may give us some computa- 
tional power using techniques of derivation or integration. 

The paper is organized as follows. Section II shows how 
a GF can find the number of distinct ways of assigning TV 
identical bosons or fermions to to different states. Section 
III discusses the GF approach to find the degeneracies of 
energy levels of a harmonic oscillator and a confined free 
particle both in higher dimensions. In our discussion, we 
follow the notation [t n ] F(t) used in Refi to mean the 
coefficient of t n term in the power series expansion of 
F(t). For example, we have [t n ] e* = 1/n! (n > 0). 



II. FINDING TOTAL WAYS OF 
DISTRIBUTING IDENTICAL PARTICLES 



The problem of how many ways there are to distribute 
N identical particles among to different states lies at the 
heart of quantum statistics. It determines the specific 
distribution function that identical particles obey, i.e. 
Bose-Einstein distribution for bosons and Fermi-Dirac 
distribution for fermions. Textbooks generally tackle the 
problem by a combinatorial method^. Here we show how 
to solve it by a GF. 

A. Identical Fermions 



Pauli's exclusion principle has that a quantum state 
can accommodate one fermion at most. This makes the 
counting problem to be much easier than that of bosons. 
The number is simply the binomial coefficient to choose 
N: 



N 



N\ (m-N)l' 



(2) 



We can do it in another way. The occupation number 
for fermions is and 1 only. We thus associate each state 
with a GF F(t) = t° + 1 1 = 1 + 1. It is obvious that the 
coefficient of term t N in the expansion of F m (£) is exactly 
the total ways that we can pick out N candidates from 
to states: (N < to) 
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[t N ]F m (t) = [t iV ](l + i) 



Ah 



1 d N 

— [t°] m(m - 1) • ■ • (to - N + 1)(1 + t) 

m ' 
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M (to — TV) ! 
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(3) 



The GF approach recovers the binomial coefficient to choose N. Furthermore, the approach can be easily 
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adapted to the boson case. 

B. Identical Bosons 

There can be an arbitrary number of bosons in one 
quantum state in principle, i.e. the occupation number 
can be 0,1,2,- • • . It is not easy to figure out a combinato- 
rial method to find the total ways in this case. With the 
help of GF, we now associate each quantum state with 
B(t) =t° + t l + t 2 H = 1 + t + t 2 H . Notice that 



B(t) has an explicit form when t is within its region of 
convergence: 

B(f) = -J_, |t|<l. (4) 



Just as what we have done for fermions, we multiply m 
copies of B(t) together. As a result, the number of total 
ways to distribute N bosons among m different states is 
just the coefficient of term t in the expansion of B m (t): 



[t N ]B m (t) 



(i-ty 



i 



1 



[t°] (-m)(-m -!)••• (-m - N + 1) 



N\ r 1 dt N (1 - t) 
1 

jv7 

(JV + m-l)! 
N\ (m- 1)! 
N + m— 1 
TV 



(-1) 



(1-t) 



m+N 



(5) 



This is exactly the same as what is derived from a com- 
binatorial method^, but we are using techniques from 
calculus for continuous variables. The GF approach also 
has an advantage when the occupation numbers are dis- 
tributed in some random way, e.g., 1,2,5,- ■ • . 

III. FINDING DEGENERACIES OF ENERGY 
LEVELS 

Energy levels are usually quantized in a quantum sys- 
tem. Furthermore, degeneracies often occur in a system 
with independent degrees of freedom. In this section we 
employ the GF to count the degeneracies. 

Problem of this kind can be formulated as follows. Sup- 
pose the system Hamiltonian H can be written as 

H = Hi + H u + H m + ■ ■ ■ , (6) 

where the subscripts I, II, III,- • ■ represent independent 
degrees of freedom. Suppose further that the en- 
ergy spectrum of Ha is {Ea {n),n} with degeneracies 
{d,A(n), n}, A = I, II, III, • • • . Since the degrees of free- 
dom are independent, the Schrodinger equation of the 
system is separable. Therefore, the energy spectrum of 
the system Hamiltonian H consists of all possible sums 
of {Ea(ti), n, A}. The degeneracy of a certain energy 
level E is then given by the total number of m-tuples 6 
(ni, ri2, • ■ • , n m ) that satisfy 



with m being the total number of degrees of freedom. 
Now we associate each degree of freedom A with a GF 

G A (t) = J2 d ±(n)t EA W (8) 

n 

and multiply them together over all A. As a result, we 
obtain a polynomial 7 of t, whose coefficients give exactly 
the degeneracies we want. To be specific, the degeneracy 
of energy level E is given by 

d(E) = [t E }l[GA(t). (9) 

A 

A. Harmonic Oscillator 

A harmonic oscillator (HO) is not only exactly soluble 
in both classical and quantum physics, but also of great 
physical relevance. For example, the problem of finding 
the Landau levels of an electron in a magnetic field re- 
duces to finding energy levels of an equivalent HO. A ID 
HO has a simple structure of energy spectrum, whose 
levels are non-degenerate and equally spaced: 

E 1D {n) = fn+^j fad, n>0. (10) 

Consider now a 2D homogeneous HO with potential 



E = E^m) + E u (n 2 ) + E IU (n 3 ) + ■ • • , (7) 



V(x,y) = -muj 2 (x 2 + y 2 ) 



(11) 
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Its energy spectrum can be solved in a planar polar co- 
ordinate system, noting that V(x, y) = V(p) = \ra\J 1 'p 2 ''. 
A much simpler way is to treat the 2D system as two 
uncoupled ID HOs. Consequently, the energy levels are 
determined by a couple of integers (n x , n y ) as 



E2F>{n) = (n x + n y + 1) hu, n > 0, 



(12) 



where n = n x + n y . For a given integer n, there cor- 
respond Ti + l different couples of quantum numbers 
{n x ,n y ), namely, (0, n), (l,n — 1), (n, 0). The de- 
generacy of energy level E2T>(n) is therefore 



1. 



(13) 



Now we apply the GF approach to find the degeneracy. 
Subtracting the zero point energy ^ hui from the spectrum 
of each ID HO and setting Hu> = 1, we obtain £ , 2d(w) = 
n x + n y , n x , n y = 0, 1, 2, • • ■ . Then we associate each ID 
HO with a GF 



(14) 



n>0 



which reduces to 1/(1 — t) for \t\ < 1. The degeneracy 
d2D(n) of energy level E2T>(n) is therefore given by 

d 2D (n) = [t n }H 2 D (t) 



= n 



i 



n + 2-l 

n 

71+ I, 



(15) 



using result from section II- B. 

We can go further to find the degeneracies of energy 
levels for a 3D homogeneous HO whose potential is 



1 



V(r) = —mLd 2 r 2 



1 



—mu> 
2 



(a; 2 + y 2 + z 2 ) 
It is straightforward to write down that 

E 3 T>(n) = 



(16) 



3 
2 

and the corresponding degeneracy 
d m (n) = [t n }H 3 1B (t) 

= n 



+ - ) huj, n>0, 



1 



(17) 



a - 1) 3 

n + 3-1 

n 



-(n + l)(n + 2). 



(18) 



The sequence of d 3 n(n) goes like 1,3,6,10,15,21,- • • , which 
is consistent with the results derived by solving in spher- 
ical coordinates. 



Energy intervals (Eq\ 
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170001-180000 


12 


7845 


3273810 


180001-190000 


11 


7842 


3366202 


190001-200000 


12 


7844 


3456580 


Total States 


447 


156634 


46596729 


Rough Value a 


4.47 x 10 2 


1.57 x 10 5 


4.68 x 10 7 



"Rough value of total states within energy E estimated by the 
quasi-continuous approximation. For ID, N = E 1 / 2 . For 2D, 
N = (tt/4)£. For 3D, N = {tt/6)E 3 / 2 . 

TABLE I: List of number of states in 20 lOOOOSo-long en- 
ergy intervals for a ID, 2D and 3D confined free particle, 
respectively. -Bo = 7r 2 fi 2 /2mL 2 . 



It is not difficult to generalize above to the case of in- 
homogeneous HOs, e.g. V{x,y) — \m{uj 2 x x 2 +u 2 y 2 ). As 
long as the ratio lo x /uj v is some rational number, we can 
always find appropriate GFs to calculate the degeneracies 
of the system energy levels. 



B. Confined Free Particle 

Consider first a ID free particle confined to some region 
of length L. It has non-degenerate energy levels^ as 



n 2 h 2 
2mL 2 



n , n > 1. 



(19) 



In the following, we choose L such that Em(n) = n 2 for 
simplicity. Accordingly, the levels Em(n) have a GF as 



F m (t) = J2 tn2 =t + t 4 + t 9 + 



(20) 



n>l 



Fxo(t) can be used to find the degeneracies of energy 
levels of a free particle confined in a square of L x L or 
even a cube of L x L x L. Finding the degeneracies in each 
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FIG. 1: List plots (dots) of Table U] which approximate the 
density of states versus energy curve (solids) n(E) = dN/dE. 
nm(E) oc E~ 1/2 , n 2 v(E) = const and n 3 n(E) oc E 1/2 . 



case (ID, 2D and 3D) is a basic problem in quantum me- 
chanics. In solid state physics-it is also concerned with 
the density of states per unit energy n(E) = dN/dE. 



In two dimensions, the density of states «2d (E) is inde- 
pendent of the energy E, whereas nm(E) oc E~ x l 2 in ID 
and n3n(E) oc E 1 ! 2 in 3D. Textbooks 9 usually derive the 
density of states by using a quasi-continuous approxima- 
tion of the distribution of energy levels in fc-space. 

Below we use Fn)(t) to find the level degeneracies in 
2D and 3D. Following Eq. (9), the degeneracies are 



d 2 B(E) 
d 3 u(E) 



E }F? D (t). 



(21) 
(22) 



Unfortunately, series like F\o(t) does not have a known 
explicit sum. Thanks to the advanced math softwares 
like Mathematica, we can use a computer to perform 
the evaluation. Mathematica's "CoefficientList" function 
can extract the coefficients in the expansion of Eqs. (20- 
22). Then we are able to count the total number of states 
whose energies fall into the interval [Ex , E2] ■ TableUgives 
a detailed list of number of states within certain energy 
intervals for a ID, 2D and 3D free particle, respectively. 
The list plots of the number of states are further pre- 
sented in Fig. [TJ There are totally 20 energy intervals 
and each interval has a length of 10000. Since the in- 
terval length is far more than the minimum separation 
between energy levels, the quasi-continuous approxima- 
tion is valid and the plots will reflect approximately the 
dependence of density of states on the energy. A com- 
parison between the discrete number of states and the 
continuous density of states is made in Fig. Q] (depicted 
as dots and solid lines, respectively). It can be seen that 
the exact numbers of states derived by a GF approach 
comply well with nm(E) oc E^ 1 ^ 2 , ri2T>(E) = const and 
n 3 n(E) oc E x l 2 . 
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